Quantum Mechanics of Spin Transfer in Ferromagnetic Multilayers 
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We use a quantum mechanical treatment of a ballistic spin current to describe novel aspects 
of spin transfer to a ferromagnetic multilayer. We demonstrate quantum phenomena from spin 
transmission resonance (STR) to magnetoelectric spin echo (MESE), depending on the coupling 
between the magnetic moments in the ferromagnetic thin films. Our calculation reveals new channels 
through which the zero spin transfer occurs in multilayers: the STR and MESE. We also illustrate 
that counter-intuitively, a negative spin torque can act initially on the second moment in a bilayer 
system. 

PACS numbers: 75.70.Ak,72.25.-b,85.75.-d 
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The phenomenon of spin transfer llUjI has resulted in 
a recent wave of progress 3, 4 0, |^ 0, i8, 9, 10, 11] to 
study interactions between spin-polarized electrons and 
a magnetic moment in a ferromagnetic film. However, 
the physics of spin transfer still has not been fully appre- 
ciated quantum mechanically. The underlying principle 
for this phenomenon is angular momentum conservation. 
Since the incoming electrons are spin-filtered by the mag- 
netic moment, as a reaction to the filtering effect the mag- 
netic moment becomes tilted to align along the direction 
of the incoming spin. This alignment is also known as the 
spin torque because of the associated moment change. 
However, spin is a quantum mechanical concept while 
torque is a classical quantity. These disparate views are 
reconciled since spin is transferred quantum mechanically 
from the incoming electrons to the magnetic moment, but 
the effect on the moment appears in a classical manner. 
In this paper, we demonstrate some quantum mechanical 
aspects of the spin transfer in a ferromagnetic multilayer 
by adopting an adiabatic approximation to describe the 
dual electron spin/ferromagnetic moment system. 

Two very interesting quantum mechanical phenomena 
in these problems are the spin transmission resonance 
fSTR) Hal and the magnetoelectric spin echo (MESE) 
effect [13. The STR is due to quantum interference 
of the right-moving (transmitted) and left-moving (re- 
flected) electron wave functions in the ferromagnetic film. 
On the other hand, the MESE is a consequence of the 
time reversal symmetry between the two magnetic mo- 
ments. The STR and MESE have some similarities as 
well as differences between them. Both phenomena oc- 
cur with effectively zero spin transfer to the moments. 
However, as we elucidate below, the physics governing 
the two phenomena is quite different. 

First, STR can occur not only in a multilayer system 
but also in a single ferromagnetic film while the MESE 
can only take place in a multilayer system. For STR, res- 
onance of the electron wave function is essential, while 
for the MESE we require an anti-symmetric configura- 
tion between the two magnetic moments all the time. In 
this way the spin current, once absorbed by the first mo- 



ment, can be generated by the second moment. That is, 
the second moment plays the role of a spin battery. Con- 
sequently, a strong interaction between the two moments 
is necessary for the MESE while STR is more or less in- 
sensitive to the moment-moment interaction. Moreover 
specific values of the kinetic energy of the incoming elec- 
trons are required for STR; this is not the case for the 
MESE. 

To account for both STR and the MESE, we consider 
a simple model Hamiltonian for two single-domain ferro- 
magnetic films with all key interactions included 
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where m is the electron mass, Jh is the coupling between 
the magnetic moment and the incoming electron spins, 
Jm is the coupling between the moments, and 70 is the 
gyromagnetic ratio. We consider the ballistic regime as 
in Refs. [ll[3. 

The magnetic moments Mi and M2 in the magnetic 
multilayers are assumed to originate from the local spins 
in the same ferromagnets with a constant magnitude Mq. 
The directions of Mi and M2 at a given time are {Oi, (pi) 
and (6*2, 02 ), respectively, and are subject to the inter- 
actions. In the Hamiltonian Eq. J^l, the second term 
transfers the electron spin to the magnetic moments via 
the spin flip process while the third term controls the 
relative orientation of Mi and M2. 

Origins of the coupling between the magnetic moments 
are the exchange interaction and a magnetic dipole inter- 
action. The effective coupling is determined by the com- 
petition between the two interactions depending on the 
distance between the magnetic moments. A simple eval- 
uation of the magnetic interaction energy shows that the 
anti-parallel configuration is more stable than the par- 
allel case when the distance between the two ferromag- 
netic films is larger than the atomic scale. For exam- 
ple, for one-dimensional uniform ferromagnets of length 
L separated by the distance d, the magnetic dipole in- 
teraction energy is approximately zL2Mq / L'^d for paral- 
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FIG. 1: (Color online) Quantum mechanical problem asso- 
ciated with the spin transfer to a ferromagnetic multilayer. 
Here rtii = Mi/Mo. Regions II and IV represent the ferro- 
magnetic films with thickness L\ and L2, respectively. Note 
that xi = Li, X2 = Li -\- d, and = Li + d + L2 where d 
is the thickness of a spacer. Regions I, III, and V are non- 
ferromagnetic. 



lei {+) and anti-parallel (— ) configurations. For two- 
dimensional uniform ferromagnets with size of L x L, 
the energy is about ±8AfQ \n{L/d)/L^. As we mentioned 
earlier, a necessary condition for the MESE is that the 
two moments remain anti-parallel for all times. Thus Jm 
must be negative with magnitude large compared to Jh- 
Otherwise, the spin torque will align easily both moments 
parallel to the initial direction of the incoming spin; the 
MESE can no longer occur. 

We show in Fig. 1 the geometry of the problem. The 
incident direction of the electrons is in the positive X di- 
rection, and the two films are placed perpendicular to this 
direction, in the YZ plane. Translational symmetry is as- 
sumed in the plane. Therefore, the problem is effectively 
one-dimcnsional as in Refs. [3, 0|; this is a reasonable 
approximation in the ballistic regime. It is also assumed 
that regions I, III, and V are non-ferromagnetic and are 
the same in nature for simplicity. When we solve for 
the dynamics of the magnetic moments we use the adi- 
abatic approximation. While this is applicable in most 
cases, it should be used with caution for the MESE be- 
cause the magnetic moments will oscillate rapidly if the 
moment-moment coupling is large. The applicability of 
the adiabatic approximation requires the time scale for 
the electrons to be much shorter than that of the mo- 
ments. 

In principle, one can use any basis set to represent the 
wave functions in each region. For example, one could 
use the spin- up/down state |±) in the lab frame for the 
basis. However, using the right basis makes the calcu- 



lations much easier. We introduce eigenstates |xct) and 
l^cr) of the interactions 2Jhs ■ Mi and 2Jhs ■ M2 such 
that 2Jhs ■ Milxa) = ±JHMo\xa) and 2Jhs ■ Majc^) 
±J//Mo|^cr)- Using the eigenstates of the interaction, 
we can represent the wave function in each region. We 
emphasize that this procedure can be systematically gen- 
eralized to a system with more ferromagnetic films. 

The wave functions in regions I, II, III, IV, and V can 
be written as follows: 



IV'///) 
IV'/v) 



-ikx 



E T.\xa){x.\+)e 



ikx 



'ikx 



(2) 



The coefficients Ra to are determined by the bound- 
ary conditions at x = 0, Li, Li -|- d, and Li + d + L2. In 
general, the coefficients depend on the directions of the 
moments in the magnetic multilayer. If different bases 
are used to expand the wave functions of each region, the 
coefficients will be changed. However, the wave functions 
remain the same. Since expressions of the coefficients are 
too long, we do not show them here. Instead, we explain 
important properties associated with the coefficients and 
the wave functions, below, for the extreme cases: i) the 
STR with no dipole interaction and ii) the MESE with a 
strong dipole interaction. 

Under the STR condition, the ferromagnetic film be- 
comes transparent to the incoming spins [l^ . One 
set of conditions that is required is that Li = 2nLoo 
(71 = 1,2,---) for an energy ratio 77 — 5/4. Here 
Loo — T^l V^mJiiMo is a typical length scale of the ferro- 
magnetic film and 77 is the ratio of the incoming electron 
energy k"^ /2m to the interaction energy JhMq. If we 
have the STR condition only for the first film in region 
II, = so that the spin state of the forward-moving 
wave in region III is the same as that of the incident 
wave. But, unlike the single layer case, ^ Q because 
the wave function in region III is partially reflected at 
X2 — Li + d and the reflected wave can then pass through 
the first film in reverse (the same STR condition applies). 
We point out that even in this case Ra ^ R'„ because the 
magnetic moments are not necessarily parallel to each an- 
other. Suppose STR takes place only in region IV where 
the second film resides. Then clearly R'^ = so that no 
wave is refiected at X2- Nevertheless, ^ because 
the forward- moving wave in region HI has both |-|-) and 
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|— ) components. Each spin component passes through 
the second film freely due to the STR condition in region 
IV. This means that the film is transparent to two differ- 
ently spin-polarized electron beams under the same STR 
condition. When we have the STR condition for both 
films, T| = Ti as well as Tj = Tj and = K = a,s 
expected. 

The other interesting limit results in the MESE. In 
this case the coupling between moments is very large, 
and antifcrromagnetic. We find Tj = T| so that the net 
spin transfer is zero because the spin state in region V 
is the same as that of the incoming electrons. It is in 
this sense that the MESE is similar to the STR for both 
films. However, if the two films also satisfy the STR 
condition, then the moments no longer rotate in response 
to the spin current, so the result is more properly thought 
of as an STR, occurring twice, once in each film. In 
all cases we have discussed so far the flux is conserved 
in each region while the spin current is not because the 
spin torque acting on the magnetic moment in region 
II(IV) is equivalent to the spin current difference between 
region I and III (and between III and V). The net spin 
torque acting on the multilayer system is the spin current 
difference between region I and V. Later, we will describe 
the dynamics of the magnetic moments in term of the 
spin current. 

The equations of motion of the two magnetic moments 
Ml and M2 are 



Fig. 2 ( Kim and Marsiglio ) 



dt 
dM2 

dt 



2JhMi X (V'j/lslV//) + JmMi X M2 
2JhM2 X (V'/vlslV/v) + JmM2 X Ml . (3) 



Note that different wave functions are used to evalu- 
ate the spin expectation values for the equations of the 
moments. We use dimensionless units as in Ref. [l^ . 
where = Mj/Aio and the time is r = jot, where 
jo = NeLoo/imSiocai) is the one-dimensional current 
with the number of incoming electrons Ne per unit length 
and Siocai = Mq/jq. In these units, the coupling con- 
stant due to the magnetic dipole interaction becomes 
a = {Siocai /'2NeLoo)JM / Jh ■ Now the equations of mo- 
tion for mi become 
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- Pimi^rriiy - ^irriix + a (nii x nij)^ 

Pi (1 - ml^) + a (riij X m.j) ^ , (4) 



where i,j = 1,2 (i 7^ j) 



7^ = i y Re [QC,t] 
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FIG. 2: (Color online) The time evolution of the magnetic 
moments (solid curves) and the spin torque (dashed curves) 
on the moments for rj = 5/4, Li = 2.6I/oo, L2 — 2.4Loo, and 
d = Loo ■ No magnetic dipole interaction is included {a — 0). 
The initial value of the moments is mi — — m2 = (0,1,0). 
The spin torque on mi^ is Ni^x ~ Qi.zx — Qiii,zx and the 
torque on m2z is N2zx = Qiii.zx - Qv.zx 



with Cia 



B'^e^^^'^. The integration range for pi and ji is given 
by the thickness of the corresponding ferromagnetic film. 

As we mentioned earlier, another way to understand 
the dynamics is using the spin current tensor; 
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where s = 1/2 is the electron spin. The spin current in a 
particular region can be calculated using the wave func- 
tion in the region. The time evolution of the magnetic 
moment is governed by the spin torque which is the net 
spin flux absorbed by the moment [9|. For example, the 
time evolution of m\z and TO2z is determined by N\zx 
and N2zx, respectively, where N\zx = Qi,zx - Qiii,zx 

and N2ZX = Qiii.zx — Qv,zx- For a wave function 
/ u \ 
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the spin current be- 
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{Re [0^02] - Re [6^62]} 
{Im [ala2] - Im [foi&2]} 



{(Kp-|a2n-(|6ip-|62p)} . (6) 



The components ai and bi are given by the wave func- 
tions in the regions I, III, and V; namely, ai (bi) is the 
forward- moving (backward-moving) component, respec- 
tively. For example, in region III, ai = X^o- ^o-l (+IXo-) 
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«2 = Y.M-\Xa){xA+), h = and 

b2=EaK{-\^.){U+)- 

We solved Eq. for a typical case with 77 — 5/4, 
Li/Loo = 2.6, L2/L00 = 2.4, and d/Loo = 1- The mag- 
netic dipole interaction is set to zero (a = 0) but a small 
interaction (|a| < 1) gives a similar result. The initial 
direction of mi = (0, 1, 0) and 1712 = (0,-1, 0). In Fig. 2, 
we plot miz and as the solid curves with the corre- 
sponding labels. Based on the single layer analogy, one 
could expect that mi and m.2 end up aligning with the 
electron spin current along the Z direction. The behav- 
ior of mi does follow this naive expectation except for 
a dip near t = 22, However, the initial dynamics of m2 
is unexpected; the magnetic moment initially acquires a 
moment in the negative Z direction, anti-parallel to the 
incoming spin current. 

We also plot the spin torque acting on the moments as 
dashed curves with labels. The behavior of the torques 
explains the dynamics of the magnetic moments. Mostly 
a positive spin torque is applied to mi so that it aligns 
along the Z direction more or less monotonically (there 
is a small dip in the torque which corresponds to the dip 
in the ruiz curve mentioned above). On the other hand 
m2 moves below the XY plane due to a negative torque 
up to r = 8. Then a positive torque begins acting on m2 
until the magnetic moment aligns to the Z direction. 

A second case is one in which a strong magnetic 
dipole interaction exists between the two films: MESE. 
In Fig. 3, we plot the numerical results for a = —500, 
7] = 15, Li/Lqo — L2/L00 ~ 2, and d/Loo = 1- We use 
a large energy for the incoming electron to ensure the 
applicability of the adiabatic approximation. The solid 
oscillatory curves with labels are the time evolution of the 
magnetic moments. The dashed line is the spin current 
in region I and V while the dashed (oscillatory) curve is 
the spin current in region III. For the MESE, the spin 
torque is not completely equivalent to the spin current 
difference because the magnetic dipole interaction also 
contributes to the torque. Nevertheless, the spin current 
still provides useful information for the dynamics of the 
moments. As shown in Fig. 3, Qi^zx — Qiii,zx is posi- 
tive. This indicates that the spin current is absorbed by 
mi. On the other hand, Qm^zx — Qv,zx is negative but 
\Qiii,zx - Qv,zx\ = \Qi,zx - Qiii,zx\- Consequently, the 
exactly same value of the spin current is generated by m2 
because of the time reversal symmetry between the two 
moments. This is why effectively zero spin is transferred 
for the MESE. 

As discussed in Ref. ^3 j the thickness of the two films 
should be the same. Brataas et al. estimated the effect 
of the thickness mismatch. We quantify this effect by 
calculating {Qi^zx — Qv,zx)IQi,zx for given parameters. 
For the same parameters except for 77 = 10, we found 
that a 2% mismatch {L2/L00 = 2.04) gives 0.3% of the 
spin transfer on the time scale of Fig. 3 (r = 20) while 
a 5% mismatch yields 1.7% of the spin transfer on this 
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FIG. 3: (Color online) The time evolution of miz and m2z 
(solid curves) for MESE with a = -500, r; = 15, Li = 2Loo, 
Z/2 = 21/00, and d — Loo. The dashed line is the spin current 
in region I and V while the dashed curve is the spin current in 
region III. The initial value of the moments is mi = — m2 = 
(1/^,1/73,1/^3). 



time scale. A longer time scale of course leads to more 
spin transfer. We also found that the spin transfer is in- 
sensitive to a (for large values), and that if 77 is increased, 
the spin transfer is decreased for the same time scale. 

In summary we have studied quantum mechanical as- 
pects of the spin transfer to ferromagnetic bilayers. Our 
formulation is readily generalized to multilayers. The 
physics of the spin transfer in the multilayers is gener- 
ically different from that of a single layer. We demon- 
strated novel quantum mechanical features which can oc- 
cur only in the multilayer system such as negative spin 
torque and magnetoelectric spin echo (MESE). The spin 
transmission resonance (STR) in a multilayer is also il- 
lustrated. Our calculation reveals new channels through 
which the zero spin transfer occurs in multilayers: the 
STR and MESE. In spite of irregularities in real mate- 
rials, we expect to see some of the qualitative features 
described here. 

This work was supported in part by the Natural 
Sciences and Engineering Research Council of Canada 
(NSERC), by ICORE (Alberta), and by the Canadian 
Institute for Advanced Research (CIAR). 
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